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Evidence continues to grow in the MiniBooNE (MB) antineutrino mode supporting a low-energy 
excess compatible with the MB neutrino mode and possibly also confirming the results of the LSND 
experiment. At least one sterile neutrino is required to explain the anomalies consistent with the 
observations of other experiments. At the same time, there is a strong tension between the positive 
signals of LSND and MB and the null results of v a and disappearance experiments. We explore a 
scenario, first proposed in jl,, where the presence of an additional heavy sterile neutrino (with mass 
well above an eV) can alleviate tension between LSND, MB and the null results of disappearance 
experiments. We compare and contrast this 3+1 + 1 scenario with the more standard 3+1 scenario 
and carry out global fits to all oscillation data including new 2011 MB v data. We find that the 
tension can be somewhat alleviated and that a phenomenologically viable window for the heavy 
neutrino, consistent with rare decays and BBN constraints, can be found if the fifth neutrino has a 
mass of order 0.3 — 10 GeV. We also find, however, that the 2011 MB v data exacerbates the tension 
with null experiments in both the 3+1 and 3+1+1 models when the lowest energy bins are included, 
resulting in little improvement in the global fit. We also discuss the implications of an additional 
neutrino for the reactor and gallium anomalies, and show that an oscillation explanation of the 
anomalies is disfavored by cosmological considerations, direct searches, and precision electroweak 
tests. 
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I. INTRODUCTION 

Neutrino masses imply the presence of new states that can generate neutrino mass terms consistent with 
the standard model (SM) SU(2) gauge symmetry. Since the observed neutrino mass splittings are tiny, the 
standard way to implement the new states is to decouple them by giving them large masses. At energies 
relevant for neutrino experiments, this gives rise to a new higher-dimension operator which generates neutrino 
masses, presumably at the scale of grand unification. Since the dynamics of the new physics is decoupled, 
however, this mechanism for neutrino mass generation can never be directly tested. 

Recent experimental hints have, on the other hand, suggested that there may be new dynamics in the 
neutrino sector at a much lower scale, leading to the possibility of probing the neutrino mass generation 
mechanism directly. The LSND [2, and MiniBooNE (MB) j3HS] experiments both have reported results 
consistent with oscillations through a new sterile neutrino mass eigenstate with a splitting that is larger 
than the splittings that control the oscillations of the SM neutrinos. The SM mass splittings are fixed 
by the observations in atmospheric and solar neutrino experiments to be O(10~ 3 eV 2 ) and O(10 _5 eV 2 ), 
respectively. By contrast, the results from the LSND and MB — > v e and — > v e searches are consistent 
with a mass-squared splitting roughly between 0.1 and leV 2 , which would require a new, heavier neutrino 
mass eigenstate. 

The simplest extension of the SM that can satisfy these requirements is a single sterile neutrino (the 
3+1 scheme). The existence of such a neutrino in the LSND and MB preferred mass region is, however, 
disfavored by global fits to the data ,6], since null searches for neutrino disappearance tightly constrain the 
mixing angles needed to produce the LSND and MB signals. Measurements of v e fluxes from nuclear reactors 
and fluxes from beam dump experiments can be combined to reject the relatively large mixing angles 
required by LSND and MB. As we show below, this statement remains true even using the new reactor flux 
predictions as inputs. Thus, a new neutrino capable of explaining the combined neutrino oscillation data 
enters a very constrained parameter space. 

In addition to these considerations, there are more complications facing the 3+1 hypothesis. Early results 
from MB [3] [5] suggested that such a 3+1 scheme might not have been compatible with the MB data alone, 
since the parameters needed to fit — > v e and — > v e appeared to be different: at high energy the MB 
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anti-neutrino mode favored oscillations and was in better agreement with the — > v e LSND data, while 
Vfj, — > v e data was consistent with a null result H The addition of a second sterile neutrino (the 3+2 scheme) 
allows for CP violation which can reconcile differences in P„ — > v e and j/„ — > ;/ e . However, this scheme 
suffers from a similar tension between the null data and the positive signals and does little to ameliorate the 
difficulties of the 3+1 scheme 0[7]. In addition, new v data [8] from MB, shown in Fig. [I] indicates that 
the apparent difference between the v and v modes is disappearing, thereby obviating one of the primary 
appeals of the 3+2 framework. 

Still, the tension between the null results and the LSND and MB data persists. In this paper we consider a 
simple scheme, proposed in [T], designed to alleviate the tension between the LSND and MB positive signals 
and the null results from reactor and short baseline experiments. This scenario requires a single light sterile 
neutrino with mass splitting in the MB and LSND range between 0.1 and 1 eV 2 and a second much heavier 
(Am 2 ^s> 1000 eV 2 ) neutrino whose oscillations are averaged over. The heavier neutrino participates directly 
or indirectly in both disappearance and appearance experiments. Because most disappearance experiments 
have their first detector relatively far from the neutrino production point the heavy neutrino has undergone 
many oscillations before reaching the detector, and the effect of the heavy neutrino is to change the flux 
of the initial flavor neutrinos. If this flux is not precisely known, as is true in many reactor experiments, 
the experiment is relatively insensitive to oscillations through the heavy neutrino. Appearance experiments, 
by contrast, look for the appearance of a new flavor in a pure initial flavor beam, so they are sensitive to 
oscillations through the heavy neutrinos. In this way, if the initial neutrino flux is not very well known in 
the disappearance experiments, new parameter space may open for the appearance experiments, giving rise 
to the possibility that the positive signals from LSND and MB are no longer in conflict with the results from 
otherwise null experiments. 

The purpose of this paper is two-fold. First, we examine the 3+1 scenario in light of the new MB v data. 
This improves the compatibility of the combined appearance data within the 3+1 framework, but we find 
that the best fit region shifts considerably to larger mixings and smaller mass splittings, which increases the 
tension with the null experiments. Second, we explore the phenomenology of, and present constraints on, 
the "3+1+1" framework of [TJ[H]. We will examine exactly how and to what extent the fifth neutrino is able 
to have an effect on the allowed parameter space of the fourth neutrino. 

The outline of this paper is as follows. We begin by establishing our notation and conventions. We then 
carry out fits for the 3+1 scenario in light of new data from MB, complete with constraints from a diverse 
set of null experiments. We then turn to discussing the parameter space for the 3+1+1 scenario with respect 
to neutrino experiments, before analyzing in detail the constraints from BBN, astrophysics and rare decays, 
which constrain the fifth neutrino to have a mass ~ 0.3 — 10 GeV. In section IV we present aspects of some 
models that explicitly realize the features of the 3+1+1 scenario, and we conclude in section V. 



II. PHENOMENOLOGY OF STERILE NEUTRINO MODELS 

We establish our notation and contrast the 3+1 framework (see e.g. [SJ [71 E] ) with the 3+1+1 scheme 
PJ [§] • We will lay out some conventions for discussing these models, leaving a more complete discussion of 
statistical methods and derivation of the oscillation formulae in the 3+1+1 scenario to the appendix. 

We aim to examine the oscillation appearance and disappearance data in depth, with a specific emphasis 
on the new conclusions to be drawn from some recently presented preliminary MB data [5] . This new data is 
in better agreement with the LSND data and prefers a sterile neutrino with lower mass and more substantial 
mixing than indicated by the earlier MB data. 



1 A low-energy excess in the f M channel was initially suspected of being a systematic effect [3] and was reported to be 
incompatible with a neutrino oscillation interpretation [I]. 
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FIG. 1: Comparison of 2010 (upper panel) and 2011 (lower panel) MB v L/E data with MB v and LSND 
v L/E data. In both panels, the MB v data is taken from [3]. In the upper panel, the MB v data is taken 

from [5], while in the lower panel the MB v data has been updated with the results of [8]. We show the 
best fit lines in the 3+1 scenario (black), the 3+1+1 scenario (green) using all data points, and the 3+1+1 
scenario (orange) dropping the three low-energy data points so that the data is in the range 
E v > 475 McV. In all plots, v lines are dotted and v lines are solid. 



A. Conventions 



We parameterize the mass mixing by 



(i) 
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where v a are the neutrino flavor eigenstates, which include the 3 left-handed (active) neutrinos of the SM 
plus any SE/(2)-singlet (sterile) neutrinos; U a i are the elements of a unitary N x N matrix that diagonalizes 
the neutrino mass matrix and causes mixing between the neutrino flavor eigenstates; and are the neutrino 
mass eigenstates with mass m, ordered by increasing mass. 

Since we are focusing on a 3+1 scheme (with a single light neutrino) and a 3+1+1 scheme (with one light 
and one heavy neutrino) we are generally interested in oscillation probabilities where all but one of the mass 



eigenstates are easily kinematically accessible. From Eq. (B3), the probability of detecting Vp in a v a beam 
is 

Pu^u, = S a p[l + 2(a-l)\U aB \\U ps \] + (l-a)\U d6 \ 2 \Up B \ !i 

- 4Es>^ ^{U* m U 0l U a3 U; j } sin 2 Xii - 4a £* =1 M{U* a5 U fj5 U aj U^} sin 2 x 5j (2) 

- 2E 8 >i>j nU* m U 0l U aj U* Pj } sin2xy - 2a £ ) =1 Q{U* 5 Up 5 U aj U*^} sm2x 5j . 



Here Xij = AmfjL/iE — 1.27 e v 2 m/MeV ' wnere L 1S the distance the neutrino has traveled and E is 
the neutrino energy. Since n 5 will be much heavier than the other neutrinos, accounting for the possibly 
suppressed production of and oscillation through n§ requires a phase space factor a that interpolates from 
(kinematically forbidden) to 1 (phase space fully accessible) as a function of the neutrino energy. For the 
short baselines and high energies of the experiments under consideration it will be a good approximation to 
take x^ ~ for i and j = 1, 2, 3, and this formula simplifies considerably. For instance, the probability for 
disappearance of flavor a is 

\-P Va ^ Va =sin 2 2^ 4 sin 2 a; 41 + 2|L/ Q5 | 2 ^l-^|C/ ct5 | 2 ^, (3) 

where we define sin 2 2Q a4 = 4|/7 Q 4| 2 (1 — |C/ Q 4| 2 — |t/ a 5| 2 ) and we assume that the characteristic oscillation 
length associated with Am^ is so short that sin 2 x^\ — > \ holds over the volume of the detector. Exper- 
iments that probe disappearance of v e are carried out at reactors and in solar neutrino searches, while 
disappearance is probed by beam dump and atmospheric neutrino experiments. 

Following pQ, the probability for v e appearance in a beam, measured by LSND and MB among others, 



simplifies. From Eq. (B4| 



with the definitions 



Pv^v^u^) = sin 2 28^ sin 2 (x ix ± ft) + k, (4) 



sin 2 2^ e = 4\U^ 4 \ 2 \U e4 \ 2 r ( . 



« - \U^\ 2 \U e4 \ 2 {(l-r) 2 + a[(l-r) 2 + 4rsin 2 p]} 
where +(— ) is for v (y) oscillations, 



(3 = itan-^i st^fel ^ l( >< 



= \u; 4 u e4 + u; 5 u e5 \ / \u; 4 u e4 \ 

sin0|[7 e5 ||C/„ B | 

and <fi = arg {jpjj^^- P is the CP-odd parameter that can account for differences in v and v oscillations. 

The 3+1 model can be recovered in the limit J7 e 5 and — ¥ 0, or r = 1 and k = j3 = 0. We emphasize 
that the sensitivity to the mixings with n 5 is such that even the limit a — > produces nontrivial oscillation 
effects. 

The 3+1+1 model is capable of opening parameter space closed by 3+1 models because of the possibility 
of CP violation and because in general we can have r > 1. In the small mixing, CP-conserving limit the 
effect of r is multiplicative because we may make the approximation sin 2 29 ^ e ~ rsin 2 26> e4 sin 2 20^/4, and 
the limits from disappearance experiments can be made compatible with larger appearance mixings if one 
has r > 1. However, we will show that because of the presence of the term that depends on |£/ e 5| 2 in Eq. ^ 
the constraints on the mixings with n 5 are almost as strong as the constraints on the mixings with n 4 . This 
forces r to be close to 1 for most of the interesting parameter space, and r is not effective in practice for 
reconciling the appearance and disappearance experiments. 
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B. Fits to Neutrino Appearance Anomalies 



Fits to the 3+1 and 3+1+1 frameworks with all relevant data are shown in Fig. [2] we display them side 
by side to enhance comparisons of the fits. In each panel we superimpose the results using the 2010 [5] and 
2011 [5] v data from MB. We use the 2009 data @] for the v mode for all fits. The best fit to the data, 
using either the 2010 or 2011 MB v data, indicates a new sterile neutrino described by a mass splitting 
Amfj - O(0.03 eV 2 ) and a mixing angle roughly of size 1, although the \ 2 is relatively shallow and is 
consistent with mass splitting Am^ ~ O(0.5 eV 2 ) and mixing angle ~ 0(3 x 10~ 3 ). These values differ 
from those we would find if we omitted the low-energy MB v and v points. As we discuss in more detail 
below, dropping these points reduces the significance of the signal so that the data are compatible with no 
oscillations at the 99% level, as noted in [8]. Because much of the significance of the fit to oscillations is 
derived from events with E v < 475 MeV, we do not omit these points in our fits. 

In principle, both appearance and disappearance oscillation experiments can bound the LSND and MB 
preferred region. We consider null appearance searches at KARMEN [T2], E776 [J3], NOMAD [HJ, CCFR 
[15] , and NuTeV [16] , and we find that the preferred region using the new MB v data is no longer in tension 
with these searches, due to the lower-mass preferred region. Although the LSND and MB oscillation results 
are not strongly constrained by the null appearance searches, the mixing angle probed by the appearance 
experiments can be tightly constrained by combining the results of v e and disappearance experiments. The 
disappearance experiments independently constrain sin 2 2# e 4 and sin 2 2# M 4, and in a 3+1 scenario with small 
mixing angles we can approximate sin 2 29^ ~ sin 2 2# e4 sin 2 28^/4, so we obtain limits on the LSND and 
MB parameter space by combining the two sets of constraints. Details of how these constraints are combined 
are given in the appendix. The v e disappearance constraints include short-baseline reactor experiments with 
new reactor flux predictions [17j^] as well as constraints from the ratio of flux observed in the Bugey 40 m 
and 15 m detectors [T5]. Disappearance of is constrained by CDHS [TU] and CCFR [2U] at high mass. 
We also take into account mass-independent unitarity constraints arising from the maximal measurement of 
the atmospheric iy^) |21j disappearance mixing angles made by the Super-Kamiokande experiment. In the 
appendix we show that this leads to + |t/ p5 | 2 < 0.0175 (0.0274) at 90% (99%) confidence. 

The best fit oscillation statistics for the 3+1 scheme are given in Table |lj where we show the x 2 nin values 
for the disappearance and appearance data sets individually as well as the Xmin f° r tne global data set. The 
value of the Xmin/DCT for the global fit does not indicate a bad fit to the data (as noted in, e.g., [10]) but 
the Xpg va l ue for the different data sets is very high, which indicates that the data as a whole are not 
compatible. This is reflected in Fig. [2] which shows the "disappearance" curve ruling out the "LSND & 
MB" allowed region. The new data play an important role in shifting the preferred region and increasing the 
tension between appearance and disappearance: we find that the appearance data on its own is marginally 
more self-consistent when incorporating the 2011 MB v data {x mm LSND+MB2011 = 24.0) instead of the 2010 



2010 Data 

Xmin bins 

Disappearance 25.4 49 
Appearance 0.20 5 
LSND + MB 32.1 30 
Everything 75.4 84 

Xpg = (Ex 2 ) min -Ex min = 17.7 
p-value = 5.02 x 10 -4 (3.48 a) 



2011 Data 

Xmin bins 

Disappearance 25.4 49 
Appearance 0.20 5 
LSND + MB 24.0 30 
Everything 72.9 84 

Xpg = (E X 2 ) min - E Xmin = 23.3 
p-value = 3.44 x 10" 5 (4.14 a) 



TABLE I: Fits to the 3+1 framework using 2010 and 2011 v data. With the new MB data, the appearance 
and disappearance experiments disagree at more than the 4er level. 



2 The new reactor flux has been reported to reflect oscillations of a sterile neutrino, but we find that it is not consistent with 
our preferred region, and we use the reactor data as a constraint. We discuss a possible resolution to this anomaly below. 
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FIG. 2: Fits in the 3+1 (left) and 3+1+1 (right) neutrino models. We also contrast the allowed regions 
using the 2010 (light-orange) and the 2011 (dark-orange) MB v data. In both panels we show the 
appearance allowed region at 99% as well as the appearance null result and disappearance null result 
exclusion curves at 99%. There is significant tension with the disappearance experiments and oscillations 
reported by LSND and MB for both the 3+1 and 3+1+1 scenarios with the 2011 data. 



data (Xmin LSND+MB2010 = 32.1), while the parameter goodness of fit becomes slightly worse (xi>G 2010 
and Xpg 2011 = 23.3). This is a result of a more significant departure from the null oscillation hypothesis at 
large L/E in the 2011 data, which is compatible with the excess at low energy found in the MB v data, as 
can be seen in Fig.JT] Due to the increased power at large L/E, our global appearance region is at somewhat 
lower mass and higher mixing than shown by previous global fits {e.g., |i6j). We conclude that the tension 
between the positive signals and the null searches indicates that a single sterile neutrino is very unlikely to 
explain the entirety of the collected data. 

In the right panel of Fig. [2] we show the results of a similar analysis performed in the 3+1+1 framework, 
where the fit region differs from the the 3+1 case because of the CP-odd phase /3 and the multiplicative factor 
r. The parameter r represents a potentially significant handle in the 3+1+1 framework, since it can give a 
multiplicative enhancement of the appearance angle compared to the disappearance angles. For small /3, r 
effectively measures the magnitude of the mixings with n 5 , and to obtain the desired enhancement over the 
3+1 scheme we need the mixings with 115, and thus the value of r, to be greater than 1. However, as shown 
in the left panel of Fig. [3j we find that r is bounded by the null experiments to be very close to 1 for most 
of the values of |/7 e 4| 2 and (t^l 2 favored by LSND and MB. This is a consequence of the "zero-distance" 
effect which allows for the oscillation of neutrino flavors at arbitrarily low distances. The zero-distance 
effect can manifest itself in two ways in the experiments in consideration. First, disappearance experiments 
bound the sum of the mixing angles due to this effect, as in Eq. (|3), which forces either |C/ Q 5| 2 or | t/ CK 4 1 2 to 
be small. Thus, disappearance experiments constrain r to be very c ose to 1 for large |[/ Q 4| 2 , as is true in the 
appearance preferred region. This effect is shown in the left panel of Fig. [3j where r is seen to be essentially 
compatible with 1 in the entire appearance preferred region. The second way the zero-distance effect would 
be visible is as a positive offset in the appearance probability. However, the appearance experiments exhibit 
transition probabilities roughly of order 0.5%, which allows them to place their own firm upper bound on 
|C^e5| 2 |^5| 2 , as is visible in the right panel of Fig. [3} In other words, we find that r is negligibly effective in 
reconciling the appearance and disappearance data sets. 

The other potential advantage of the 3+1+1 framework is the possibility of CP violation, but we also find 
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FIG. 3: Constraints on r as a function of |C/ e 4| 2 |i7 At 4| 2 for 0.04eV 2 < Am 2 ^ < 6eV 2 . We see that r is close 
to 1 in the appearance preferred region, and has limited ability to reduce the tension with null experiments. 



that this is not very effective in reducing the tension with the null experiments. As shown in the left panel of 
Fig. [4j the x 2 nas a pronounced preference for a small nonzero value of f3, and the sharpness of this feature 
means that the extra parameter freedom is largely unimportant in defining our preferred region. When we 
drop the MB v and v data points below 475 MeV, as advocated in the initial MB data release [3J, we find 
that CP violation has a much more significant impact on the fit. This is because the % 2 is substantially 
flatter as a function of /3 and exhibits two rather broad and nearly degenerate minima, as shown in the 
right panel of Fig. [4j This in turn is a result of the mostly flat spectrum of the v and v data points in the 
region E v > 475 MeV, which can be made compatible at the 1 a level for a wide range of low-mixing-angle 
oscillations when CP violation is allowed. The ultimate effect of the shallower x 2 is to reduce the preference 
for a particular mass or mixing, which opens up a wider range of parameter space and broadens the preferred 
region. When we perform fits analogous to those in Fig. [2] for the data with E v > 475 MeV, we find that the 
significance of the signal drops so drastically that the remaining data are consistent with no oscillations at 
the 99% level, as noted in [5]. We show the CP- violating best fits to the E v > 475 MeV data alongside the 
best fits to the E v > 200 MeV data in Fig. [l] 

Finally, we give the best fit statistics for the 3+1+1 model in Table [11] taking E v > 200 MeV as usual. As 
in the 3+1 case, the 2011 MB v data provides slightly more agreement in the combined appearance data than 
the 2010 data. Again, the global fit to the data gives an acceptable x 2 /DOF, but the PG test underscores the 
point that the data sets are incompatible. The p-value for the Xpg m the 3+1+1 case is slightly lower than 
in the 3+1 model for both the 2010 and 2011 MB v data. This suggests some improvement in agreement, 
but with the new data the tension remains at the 4cr level in both the 3+1 and 3+1+1 cases. 



III. PHENOMENOLOGY OF THE HEAVY NEUTRINO n 5 IN THE 3+1+1 SCHEME 



We discuss the parameter space of interest for the heavy neutrino, 715, that appears in the 3+1+1 frame- 
work. If rc.5 is a Majorana neutrino the neutrinoless double beta decay constraints are extremely restrictive 
[2"4"j . so we will take n.5 to be a Dirac state. 
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FIG. 4: x 2 as a function of the CP-odd parameter /3, with (left) and without (right) the MB v and v data 
for 200 MeV < E v < 475 MeV. We show fits utilizing both the 2010 and 2011 MB v data, and we show the 
90% or 99% allowed value from the A\ 2 test that we use. 



2010 Data 
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Disappearance 24.6 

Appearance 0.20 

LSND + MB 28.3 

Everything 74.4 



bins 

5 
30 
84 



xpg = (E x 2 ) min - E xLin = 2i.; 

p-value = 1.6 x 10~ 3 (3.16 a) 



2011 Data 



Disappearance 
Appearance 
LSND + MB 
Everything 



Xmin 

0.20 
19.4 
73.2 



Xpg = (Ex , 1u 
p-value = 6.0 



min 

X 



bins 

5 
30 
84 



- E Xmin = 29.0 

10~ 5 (4.00 a) 



TABLE II: Results of fits to the 3+1+1 framework using 2010 and 2011 MB v data. With the new data, 
the appearance and disappearance data sets still disagree at about the 4cr level, with only slight 

improvement over the 3+1 case. 



We begin by showing how some experimental anomalies recently reported at low significance could plausibly 
be explained by the existence of this heavy mass eigenstate. Then we proceed to place constraints on its 
parameter space. Tension with BBN constraints forces us to the regime where m 5 > 1 MeV. Above this 
region constraints from SN1987A, pion and kaon decays, beam dump experiments, and nonobservation of 
/i — > e7 enter. We find that if one wishes to use n§ to explain experimental anomalies, only a small window 
with TO5 ~ 0(1 GeV) and mixing | C/e5 1 1 1 ~ O(10~ 2 ) is allowed, with the additional requirement that 715 
be stable on collider timescales or decay to non-SM final states. 



The Gallium and Reactor Anomalies 



In recent years, several experiments have reported observing anomalously low neutrino fluxes. Since 
these anomalies include many different energy and distance scales and exhibit no L/E dependence, we do 
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not include them in the fits to 774, whose mass splitting will cause visible oscillations at these experiments. 
Instead, we fit to this data with the heavy neutrino 775 whose oscillations are averaged over in all experiments. 
A somewhat more detailed discussion of the anomalies we fit with n 5 oscillations is given in the appendix. 

The most statistically significant of these position- and energy-independent anomalies is the reactor an- 
tineutrino anomaly [17J [25] (RAA) , where the global average of the observed v e flux is less than anticipated 
by a factor -Rraa = 0.943 ± 0.023. In addition, anomalously low measurements of the v e scattering gallium 
cross-section |26j . the gallium anomaly (GA), may indicate disappearance of v e . The deficit based on four 
measurements of the average v e scattering cross section from the process v e + 71 Ga — > 71 Ge + e~ is, with 
correlated errors taken to be those in [T7J, Rga = 0.86 ± 0.06. Independent measurements of the strength 
of the relevant Gamow- Teller transitions [57] supports the conclusion that this deficit might be due to av- 
eraged oscillations of a heavy neutrino. Finally, measurements of the energy dependence of the scattering 
cross section in the process 24 + 12 C — > 12 N g . s . + e~ [28] are very mildly discrepant with the cross-section 
predictions |29j and might also be due to a similar reduction of flux of v e . This was originally presented as a 
constraint on the GA parameter space in [3D] , but the shallowness of the respective A% 2 's and the similarity 
in the parameter space leads us to consider the possibility of reconciling this data with the RAA and GA 
data. 

Since the disappearance formulae are not sensitive to CP violation, the probability for both v e and v e 
disappearance is given in Eq. ([3]). In all of the experiments in consideration, the neutrino energy is less than 
what we will find is the allowed range for 7715, so we set a — > in all cases. It is clear from Eq. ^ that 
77.5 can reduce the v e flux by a fixed amount with no energy or distance dependence whether or not 775 is 
kinematically accessible. To extract the most conservative limits on the mixings with 775 we will set sin 2 2f? e 4 
to in these fits. We fit to all of the available data using correlation information as in the literature. We 
find a preferred value of 

\U eb \ 2 = 0.036 ±0.013. (7) 

Thus, we find that the RAA and GA may be consistently reconciled with the carbon data in the presence of 
a very heavy neutrino with averaged oscillations and a mixing angle of the magnitude indicated by Eq. ([7| . 
However, we will show below that for a very massive sterile neutrino, a mixing angle of this magnitude is 
disfavored by a combination cosmological considerations, direct searches, and precision electroweak tests. 



B. Big Bang Nucleosynthesis 

Depending on the mixing and mass, additional light sterile ncutrino(s) (with mass < 1 MeV) can be 
thcrmalized in the time leading up to big bang nucleosynthesis (BBN). The presence of additional neutrinos 
at BBN can drive a faster expansion rate, modifying the abundance of the light elements, and in particular 
of helium. Thus, detailed observations of primordial elemental abundances from BBN can constrain the 
properties of sterile neutrinos. To set bounds, we require that the total number of neutrinos at BBN is 
less than 4.4 [31], or AA„ < 1.4. With potentially two sterile neutrinos with masses below 1 MeV in the 
3+1+1 scenario, BBN constraints must be carefully checked. While constraints from BBN can be alleviated 
by the inclusion of a large lepton asymmetry (which effectively delays the time when an MSW-like coherent 
conversion can occur) |32j . this mechanism becomes ineffective for the large mass splittings of interest for 
the 3+1+1 model. We review the constraints in this section and apply them to the 3+1+1 scenario. 

We begin this discussion by reviewing the calculation for one active plus one sterile neutrino. With this 
result in hand, we will be able to easily see how the result extends to two sterile neutrinos (with widely 
separated masses) mixed with more than one active neutrino. We follow the density matrix formalism of 
[33j . Assuming that the active neutrinos are always in a fully thermalized state, the evolution equations for 
an arbitrary number of neutrinos N is 

ip = [H,p]-iTp, (8) 

where the Hamiltonian T-L a p = V a p + Y^iLi UaiU^mf /2E, and we take the production rate to be T a p = 
(r„ + Tp)/2, with Vi S = Ts = for all i. Specializing to the case of one active and one sterile neutrinos, the 
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relevant equations are 

Hxd x p ss = iH as (Pas - Psa) (9) 
Hxd x p as = -i[{Uaa - U ss ) - ijas]Pas + iUasiPaa - Ps 



J SS J ; 



with x = m/T, and m fixed to be lMeV. The effect of interactions encapsulated in 7 as is to damp away 
the coherent off-diagonal element, p as - Thus if 7 as is large we are forced into the stationary point where 
d x pas » [33], so that 

Pas = ; (Paa ~ Psa)- (10) 

Substituting this in the differential equation for p ss we obtain 

„ R _ 7a , _ , sin 2 261 

tlXa xPss - 4 \J> m Pss)( cos2& _ Vaa/dE) 2 +J 2 /ASE 2> 

where SE = Ato 2 /2£, 7a = 3a ^|^G|,r 4 p and V OQ ~ -C a G%T*p/a, with 5 „ c ~ 4, ~ 2.9, C„ e 

0.61, C„ iUr ~ 0.61. If we neglect the 7^ term in the denominator, which is valid for the non-resonance case, 
this is easily soluble analytically, since the result takes on the simple form: 



la (T = TO) 



ln(l - AN V ) « sin 2 26* / dx— c — - -, rr^, (12) 

V v > 4H{T = m) J (x 6 cos26- V aa (T = m)/6E(T = to)) 2 ' 1 ' 

Doing the integral analytically or numerically, we see that the result scales as ~ ^J8E(T — m)/V aa (T = to) ~ 
V Ato 2 . The physical meaning of this result is clear. The rate with which the sterile neutrino is populated is 
suppressed at late time because the interaction rate is dropping as 1/x 5 . At the same time, the sterile neutrino 
is most likely to be populated when the mass splitting between the active and sterile states is smallest. The 
medium dependent mass splitting, however, is also dropping with V aa . Altogether, the integral is dominated 
by when x 6 cos 28 ~ \V aa (T — m)/8E(T = m)|, so that the total result (squared) is [33 

(Am^/eV 2 ) sin 4 26 es = 3.2 x 10~ 5 ln 2 (l - AN V ) (13) 
(Am^/eV 2 )sin 4 26i M , TS = 1.7 x 10~ 5 ln 2 (l - AN V ). 

With these physical insights, it is easy to see how the results generalize to the cases with more than one 
sterile or active neutrino. First, we can see that because the sterile neutrino is populated around when 
Vaa/SE ~ 1, at any given temperature only one of the sterile neutrinos will be populated if the masses of 
the sterile neutrinos are widely separated in mass from each other and from the active neutrinos. Thus, if 
we make the assumption that m a -C TO4 <C 7715, we can decouple the fourth and fifth neutrinos from each 
other and treat them as being populated only through their interactions with the active neutrinos. 

The other complication to consider is mixing between the active neutrinos themselves. However, if the 
active neutrino mass splittings themselves are negligible in comparison to the sterile neutrino mass splittings, 
Ato 2 2 <C Ato-23 <C ATO34 <C Amis, then the SM mixing angles 812, #13, 823 can be rotated away, and as a 
result the mixing between the active neutrinos itself decouples. Thus we conclude that the constraints on 
active-sterile mixing can be decoupled accordingly, and we have 

(Am 2 4 5)1 /eV 2 ) sin 4 26» e4 , 5 = 3.2 x 10~ 5 ln 2 (l - AN U ) (14) 

(ATO 2 4 5)1 /eV 2 ) sin 4 26^4,5 = 1-7 X 10~ 5 ln 2 (l - AN V ). 

Now LSND and MB, in the standard 3+1 scenario, probe sin 2 28^ = 4|C/ e4 | 2 |?7 Al4 | 2 ~ \ sin 2 26> e 4 sin 2 26> m4 , 
which is a good approximation in the small mixing angle limit. Thus we are able to conclude that 

4(Ato 2 4 5)1 /eV 2 ) sin 2 28 ^ = V3.2 x 1.7 x 10~ 5 ln 2 (l - AN U ), (15) 

so that we learn that the sterile neutrino is thermally populated if it has mixing angles large enough to 
explain LSND plus MB in the 3+1 scenario. 
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Now these results are easily extended to the 3+1+1 scenario. Then we have constraints on r sin 2 29^ e = 
4|C/* 4 t/ A1 4||C/* 4 f7 M 4 + J7* 5 J7 Al 5| from LSND plus MB. To alleviate the constraints from the disappearance exper- 
iments, we require {U^U^] > \U* A U^4\. Since m 5 m 4 , we conclude from Eq. (15) that if j/ 4 is populated, 
then v§ is also populated at BBN temperatures, unless 7775 > 1 MeV. Since AN„ < 1.4, we thus conclude that 
most of the 7775 parameter space proposed in pQ is not consistent with the constraints from BBN, eliminating 
the entire region of parameter space with 33 eV < 7775 < lMeV. 

Above the upper end of this mass range production of heavy sterile neutrinos may be inefficient at BBN 
temperatures, and the additional neutrino n 5 may not represent a fully populated degree of freedom. Because 

the state ?7 4 is fully populated, we require the fractional population of the state n$ to satisfy ANjf 1 ^ < 0.4 as 
calculated from Eq. ( 15 1. After their production is frozen out, the remnant 775 will decay through charge- and 
neutral-current interactions at rates suppressed by the mixing parameters. If these decays proceed through 
SM channels, the decays of 775 can add entropy and ionizing energy to the thermal bath at the time of BBN, 
spoiling predictions of the relic helium abundance. These considerations allow us to rule out the range 



6.2 x 10- 10 ( 1 < \U e5 U u5 \ < 10 15 (^- ) 3 5 , (16) 
VMeV/ ~ 1 M 1 ~ VMeV/ V ; 

where the upper bound is a rough fit to the numerical analysis conducted in [33] (we display the numerical 
values in Fig. [5|. These bounds extend up to 7715 ~ m n , at which point new decay channels open which have 
not been analyzed numerically. 

In the next section we consider further constraints from colliders on such heavy neutrinos. We will find 
that to satisfy the constraints, 71.5 must have exotic invisible decays which are not via its SM mixing with 
the active neutrinos. 



C. Supernova 1987A 



The duration of the observed neutrino burst from Supernova 1987A (SN1987A) constrains the mass and 
couplings of any massive sterile neutrino. If the sterile neutrinos mix strongly enough that they are produced 
but are coupled weakly enough that they arc not tightly bound to the supernova core they will allow too 
much energy to escape from the core, reducing the observed duration of the blast. 

There are both lower and upper bounds on the neutrino coupling 35J. The lower bound comes from 
requiring that 715 are efficiently produced in the interior of the supernova. If these neutrinos arc efficiently 
produced and have low enough mixing, they will free stream out of the supernova and conduct energy away 
from the core too quickly. With larger mixing angles, the neutrinos will have a short mean free path and, 
for large enough mixing, they will be trapped in the supernova. If they are trapped but their mean free 
path is larger than the supernova core they will cause anomalous cooling of the star: blackbody radiation 
will be emitted from a region larger than the supernova core, and the supernova will cool too quickly. 
Because production of neutrinos in supernovae are dominated by charge-current processes, we find that u e 
production dominates production |36j . This gives slightly weaker bounds on the mixing angle U^, and, 
because maximal mixing angles are in principle allowed by these arguments, we find that the lower bound on 
Up,5 (which is approximately 5 times weaker than the lower bound on U e §) is in fact the lower bound on the 
product {Ue^U^l. For m 5 > 0.1 MeV, where the matter effect becomes unimportant, we find that mixing 
angles 3.0 x 10~ 5 < | C/e5 C^t5 1 < 5.0 x 10~ 3 are ruled out by these energy considerations. These bounds apply 
to 715 regardless of its couplings. 

The trapping argument given above will not apply for large mixing angles if 775 decays invisibly to products 
with no SM interactions, as naturally considered in the model building section below. This is because, for the 
widths calculated below, we find that the decay length L ~ 10~ 10 m (10 MeV /m 5 ) is much shorter than the 
mean free path A m f p ~ 0.1 m / sin 2 29 m for the masses and mixings of interest. Therefore 775 will decay well 
before it is trapped, and the bounds can no longer be lifted at very large mixing angles. These exclusions are 
model-dependent because they rely on the unknown couplings of the decay products, but since it is possible 
that 775 evades the upper bounds described at large mixing angles we shade this region gray in the left panel 
of Fig. [5] 
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These bounds are also lifted for lighter sterile neutrinos since the relevant production mechanism is matter- 
enhanced flavor transitions. For m < 0(0.1 MeV) the bounds weaken and go to zero around m ~ 0(100 eV) 
[55] . so SN1987A bounds do not constrain n^. 



D. Bounds from Light Mesons 



There are a variety of searches for exotic meson decays that produce strong bounds on the mass and 
mixing of n 5 . We group these into a few categories as follows and display the collected results in Fig. [3J 

Measured it meson branching fraction: The pion branching ratio to /x and e is = |-M|^ ejy /|A1|^„ W . 
At tree level, the matrix element |A4|^._^ „ goes like 

N a 

\M\l^t a n « £ P m \ 2 (m 2 a + mj) [ml - (m a + m,) 2 ] , (17) 

i=i 

where m a is the mass of the charged lepton l a . In the SM, where to.; = for all neutrinos and Ue a i = 5i a i, 
i?7r simplifies considerably. The most current SM calculation of this quantity to two loops is ijSM,th. _ 
(1.2352 ± 0.0001) x 10" 4 [37], while the best experimental bounds give R^ = (1.230 ± 0.004) x 10" 4 [3S]. 

In the 3+1+1 framework R v will differ depending on the mass range, so the constraints are piecewise. 
They simplify at high mass, where TO5 > to^ — to m 3> m e , which is near where the SN1987A bounds stop. 
We find 



R TV 



l-\Ues\ I [gg5T m i m Z- m 5 _ < < 



oSM,th. - i-it/. a r TOk > 



(18) 



The measured ratio is i? w /i?™' th - = 0.996 + 0.003, so at 99% confidence we require that i? w /i?| M - th - < 1.004. 
Thus, the mixing angles U e $ and are bounded fairly strongly in the intermediate mass range. We do a 
scan over the full parameter space and for each value of m§ we find the maximum product of the mixing 
angles consistent with this constraint. 

Muon lifetime: For nonzero U e 5 and U^, the total charged current interactions with the muon and 
electron below the muon mass will be reduced. The muon lifetime will be increased relative to the SM 
prediction due to the non-unitarity in the neutrino mixing matrix. In practice, Fermi's constant, Gp, is 
measured most precisely from measurements of r M [39] . so one can derive constraints by comparing to an 
independent measurement of Gf- Following [40], we relate Mz, Mw and a to Gf by 

Ql = naM Z / lgl 

F j2M w {M z -M w ){l-&r) 

where Ar = 0.0362 ± 0.0005 [35] is the correction to the trcc-lcvcl relationship. The values of Mz and Mw 
used should be taken from purely kinematic measurements since other fits to M\y include the measurements 
of G F from muon decay. We take M w = 80.387 ± 0.016 GeV [41] and M z = 91.1875 ± 0.0021 GeV [42]. 
Plugging in these values, we find G' F = (1.1679 ± 0.0013) x 10~ 5 GcV~ 2 . Comparing this to the value 
extracted from measurements of t m , Gf = 1.166353(9) x 10 -5 GeV -2 [33], we find for to 5 > to m 



Gf 
G' f 



resulting in an upper limit 



(1 - |C/ e5 r)(l - \U^\ 2 ) = 0.9987+ 0.0011, (20) 



lUesU^l < 0.0021 (21) 
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at 99% CL. We mark this line as t m . 

Searches for lines in tt and K meson decays: Measurements of tt, K — > £ a n give important bounds on 
the mass and mixing of 775 with v a . These are summarized in Figs. (2-4) of |24j . Bounds are given by null 
searches for peaks in the spectra of the leptonic products of these decays. For n 5 produced by the decay of a 
heavy parent particle M of mass m« with a decay partner l a of mass m a we expect to see a monochromatic 
line in the lepton spectrum at E a = {m 2 M + tto^ — mj) /2rriM- These lines are generically not found, and 
limits on 775 mixing are based on the specifics of the given experiment. 

For the electron neutrino sector, the decay 7r — > ev 03] is strongest below m w and K — > ev [U] is strongest 
between m T and tuk- For the muon neutrino sector, the important decay is K — > fiu |45j . In the region 
m 5 > — rria, muons cannot be produced in tt decay so there are no muon bounds in that range. Thus, 
line searches of tt decays do not provide strong constraints on the product It/es^sl i n the mass range 
ms > m n — m M since experiments cannot set any bounds on in this range. 

Decays of n 5 : If n 5 is heavy and can decay to SM products, these decays will be seen in dedicated 
searches such as, e.g., the PS191 CHARM @7], and DELPHI gS] experiments. PS191 looked for the 
decay of a heavy neutrino through a variety of weak interaction channels; it is constraining from ~ 1 MeV[^] 
to 138 MeV. CHARM searched for decays 775 — > l + t~v : where I = e, fi, with constraints from 500 McV 
to 2.8 GeV. DELPHI also looked for a wide variety of 775 decays, and it provides limits from 2 GeV to 90 
GeV. We show these excluded regions as well as the limits from 775 decays in dileptonic K decays |35j . Note 
that, as pointed out previously [49], the PS191 and CHARM collaborations considered 775 decays through 
charge-current channels only. When the necessary neutral-current contributions are added |24] , the bounds 
are strengthened somewhat compared to the published results 05]. We provide bounds including both the 
charge- and neutral-current contributions. 

Non-observation of /i — > ej: For the decay \i — ► e-f, we have the standard result 



where g(rn,i) is a kinematic factor given in [24]. This is not constraining below 0(1 GeV) and by 300 GeV 
the bound asymptotes to | t/e 5 5 1 < 5.25 x 10~ 5 using the current measurement Br(/i — > ej) < 2.4 x 10 -12 
[50] - At high mass, this is the most important constraint. In particular, measurements at the Z-pole are 
weaker than /j, — ¥ &y, so we do not show these bounds on our plots. 



In Fig. ml we show bounds on the product lU^U^l, which in the CP-conserving limit is the product that 
sets the value of r in the appearance probability formula, Eq. Q. There are several model- independent 
bounds: as described above, BBN is most constraining below ~ 1 MeV; there are universally constraining 
bounds for masses 0.1 MeV < TO5 < 100 MeV from SN1987A; the NuTeV oscillation search [TB] rules out 
mixing angles |J7^5t7 e 5| > 1-3 x 10 -2 for the entire mass range; and for masses 7715 > 64 GeV, the bounds 
from /i — > e7 are most stringent. In the range 100 MeV < m$ < 64 GeV the constraints bifurcate depending 
on whether n 5 decays to charged leptons or remains invisible on collider timescales. 

Invisible decays: When n.5 remains invisible on collider timescales there are constraints from line searches, 
the pion branching fraction, and precision electroweak measurements of Gf- From m T — < 777,5 < m Tr 
the strongest bound is from the measured branching fraction of pion decays, For < m 5 < mx — 



3 The PS191 experiment did not publish limits for mixing angles above 10 4 , so we extrapolate the bounds down to 7715 = 2m e 




2 



(22) 



E. Combined Bounds on 715 



as a power law with oc 
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FIG. 5: Exclusion regions from BBN [31] (right frame) and SN1987A [35] (both frames), as well as bounds 
from the NuTeV oscillation search [16] (red dotted, right frame), R v [37j|38j (left frame), measurements of 
Tp [55HI2"] (left frame), collider and line searches [551 H3"HI9"] (both frames), and searches for // — > [50] 
(both frames). The left panel shows lines of constant values of r from 1.05 to 2.4 (for the calculation of r, 
we assume no CP violation and take |C/ e 4^4| = 0.023, as explained in the text). To avoid clutter, we avoid 
repeating the r M and NuTeV lines in both plots, although each is valid in both cases. 



searches for leptonic lines in kaon decays are constraining for both e and /i products. For m„ < m 5 , 
comparing the values of Gf from measurements of r M and the W and Z masses as described above gives 
tight constraints. We show these bounds in the left panel of Fig. [5} 

Between the K line searches and the fi — > curve, where 387 MeV < 777,5 < lOGeV, the most con- 
straining bounds on 77,5 come from the precision electroweak measurements of Gf- Although this is the 
least constrained region, we find that these measurements still disfavor large values of r. Assuming no CP 
violation and taking lE/^L/^l = 0.023, which is the smallest value of {U^U^l for which \U e ^U^\ can take 
on arbitrarily low values in the MB and LSND region, we find that r < 1.09. 

Visible decays: In addition to the SN1987A and [i — > bounds and the low mass constraints on 
BBN, we find that the direct searches at PS191, CHARM, and DELPHI are very constraining if 715 decays to 
SM particles on detector timescales, and we also find that the BBN constraints can be extended to 7715 ~ m^, 
as described above. These give the most powerful constraints from < 1 MeV to 64 GeV. Above this range, 
the \i — > e7 constraints become powerful. These bounds are in the right panel of Fig. [5] 

We see in Fig. [5] that the bounds are prohibitively strong if 715 decays to SM products. We find that 
77,5 is phenomenologically more viable provided the decays of 775 are invisible and the mass satisfies 
387 MeV < m 5 < 10 GeV. However, when restricting the range of lU^U^l from the muon lifetime, 
the LSND and MB results strongly favor r ~ 1. Furthermore, the combination of constraints from SN1987A 
and the muon lifetime restricts | C/ e s | 2 < 0.004 for 7775 > 100 MeV, which seriously constrains the parameter 
space for solving the RAA and GA data, as indicated by Eq. 

In the next section we construct models of neutrino mass that naturally allow for invisible decays. 
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IV. NEUTRINO MODELS 



The 3+1+1 scenario relies on the presence of a heavy neutrino with a substantial mixing with the light 
neutrinos. Within the standard see-saw scenario, with one active neutrino v a and one sterile neutrino v m , 
this is not possible to achieve. The mass matrix 



M 



mo 
m D M 



(23) 



connects the mixing to the mass hierarchy, so that a heavy sterile neutrino necessarily has a small mixing 
with the SM neutrino: 9 ~ mo /M, which is small for a sizable neutrino hierarchy. 

A large mass hierarchy and a large mixing can, however, be achieved for a Dirac sterile neutrino. Taking a 
single active neutrino v a and a sterile neutrino v d with Dirac partner D d , we can write a general mass matrix 
in the (v a , Vd, Vd) basis as 



M 



mi m o 
mo m 5 
m 5 



(24) 



Defining M 2 = m\ + m D and expanding to second order in the small ratio mi/m 5 we find the eigenvalues 



A 2 = —M 



($) - 



m L> ( m D + 4mg) 

mf> (m D + 4m§) 
8M5 



m( + ... 



m\ + 



(25) 



corresponding to the (unnormalized) eigenvectors 



M 



(1 + ^) 
1 



#2 = 



M 



D (-1 _ mi' 
I \ ± M . 



m 5 
M 



m 5 
M 





( 


M 2 


\ 






rn i m d 




j K 3 = 




l 






{ 


M 2 


/ 



(26) 



if i and K<i correspond to the components of the mostly sterile fifth mass eigenstate n§ , whereas K% corre- 
sponds to a mostly active light state. The mixing between n 5 and the light state is controlled by m D /M. 
This ratio need not be very small since the mass of the light neutrino is fixed independently by m\. The 
small mixing scenario is recovered in the limit mo <C m§ , which corresponds to m§ — > M, while maximal 
mixing corresponds to the limit m D — > m 5 . 

This type of scenario can be extended to encompass the fourth neutrino, as well as the needed invisible 
decays of 715 . Consider adding to the Lagrangian a term 



£«* = \4>v d v m + A'c^m- 



(27) 



Neglecting Majorana mass terms for illustration, we find that a mass matrix in the {v a , v d , v dl v m ) basis with 
the desired phenomenology is given by 



M 



( m u \ 

mo m§ m0 

m 5 

V to / 



(28) 



nr, , + m\. In 



where = Xv^. This matrix has two zero eigenvalues, with the other two set by ±^/m^ 
a hierarchy where m^ <C mo -C m 5 , the massive states are mostly Vd and their mixing with the active 
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neutrino is controlled by mu/m^. The massless states are predominantly composed of v a and v m , and their 
mixing is controlled by m^/mp. Of course, these masses should not exactly vanish, and the masses can be 
lifted from being zero by appropriately small Majorana mass terms. 

The new state <fi allows both for large Vd — v a mixing and for invisible decays of Vd (via v& — > <j>v m with 
subsequent decays (p — > v m v m ). This decay, with width ~ -^\ 2 m^^ should be compared to the decay 

derived from mixing with active states, which scales as T$m ~ ltw^i e ^z (mz) m5 ' Since m 5 ^ s m * ne GeV 
range, the SM decay channel is naturally suppressed with respect to the invisible decay. 



V. CONCLUSIONS 



We have studied fits to the LSND and MB experiments within the context of 3+1 and 3+1+1 [T] scenarios. 
Compared to the 3+1 scenario, the 3+1+1 framework posits that the presence of an additional heavy 
neutrino which is not directly probed by most disappearance experiments lifts some of the constraints of the 
null disappearance experiments. However, using the new 2011 MB v data, we find there is still significant 
tension between positive and the null results, even with the additional very heavy neutrino. We went on to 
explore the phenomenology of the massive neutrino that appears in the 3+1+1 scenario, and we found that 
a heavy mostly sterile neutrino could be consistent with a variety of cosmological and collider constraints if 
the sterile neutrino has a mass around a GeV and does not couple primarily to the SM. We also showed that 
in the face of BBN, direct search, and precision electroweak bounds, even a heavy state that decays invisibly 
might not be suitable for reconciling the anomalous measurements of v e fluxes made by gallium and reactor 
experiments. 

If the larger mixing angle required by the 2011 v data for the LSND and MB anomaly persists, other types 
of scenarios will be required in order to obtain a consistent global explanation of the neutrino oscillation 
data. One possibility, which was explored in [52], is to make use of medium dependent neutrino masses 
[53"] . In this case, Bugey (whose oscillations would mostly occur through air) would be weakened relative to 
LSND and MB (whose oscillations mostly occur through earth) , and a wide swath of parameter space would 
remain. We leave this possibility for future consideration. 
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Appendix A: Fit Details 



We summarize and expand upon the details of our analysis of experimental data included in our fits. 

Appearance Experiments: The characteristics of the experiments used in our fits are shown in Ta- 
ble "~ 



III 



LSND [2] observed the appearance of v e with energies 10 — 60 MeV in a beam of v^, consistent with 
neutrino oscillations that occur in the Am 2 ~ 0.2 — 10 eV 2 range. MB also measured P M — > D e oscillations 
[5J with energies 200 — 3000 MeV consistent with evidence for antincutrino oscillations from LSND. MB 
did not initially report evidence for oscillations of the form — > v e [2j , but an in-depth analysis published 
after the release of the initial data set supported the interpretation of a low-energy excess consistent with 
oscillations @J. We use all of the MB v data points in our analysis, including those below 375 GeV which 
were excluded in the first MB analysis. Due to the low energies of these experiments, we take a = for our 
fits. 

We also include the null results of KARMEN [12], E776 [33], NOMAD [14], CCFR [15] and NuTeV [16]. 
For each null experiment, we find a single data point-that is, the oscillation probability and error for the 
L/E value-which best matches the 90% exclusion curves given by the experiments. Due to the generally 
high energies of these experiments, we use a = 1 for all experiments except KARMEN in these fits. 
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Experiment 


mode 


# points 


Distance (m) 


E 


Am 2 (eV 2 ) 


MB 




11 x 2 


541 


200 - 3000 MeV 


£ o.i 


LSND 


Up 


8 


29.8 


10 - 60 MeV 


£ 0- 3 


KARMEN 


v„ 


1 


17.7 


1 - 50 MeV 


> 1 


E776 




1 


1000 


1 - 10 GeV 


> 1 


NOMAD 




1 


625 


> 10 - 200 GeV 


> 10 


NuTeV 




1 


1436 


> 10 - 300 GeV 


> io 2 


CCFR 




1 


1436 


> 10 - 300 GeV 


£ io 2 


TOTAL 




30 pos., 5 null 


~ 10 - 1436 


10 MeV - 600 GeV 


> o.i 



TABLE III: Energies, mixings, and mass splitting sensitivities for each appearance experiment. 



Disappearance Experiments: The v e disappearance constraints include short-baseline reactor exper- 
iments with new reactor flux predictions |17j plus constraints from the ratio of the flux observed in the 
Bugey 40 m and 15 m detectors |18) The statistics of the constraint on v e disappearance is dominated by 
the Bugey ratio. Disappearance of is constrained by CDHS [19] and CCFR [20], which we take as single 
data points, corresponding to the combined oscillation probability for the full energy range. Because both of 
these experiments search for muon neutrino oscillations between two detectors, very large mass differences 
are not restricted, since the beam is likely to be fully oscillated as it arrives at both the near and the far 
detector. Also, because of the baselines and energies of these experiments, it is a good approximation to 
ignore the probability of oscillation through n 12j 3. We take a = 1 for CCFR and CDHS, but have a = for 
the reactor experiments. 



Experiment 


mode 


# points 


Distance (m) 


E 


Am 2 (eV 2 ) 


CCFR 


v„ 


1 


714 and 1116 


40 - 200 GeV 


10 - \tf 


CDHS 


v» 


1 


130 and 885 


2-6 GeV 


10" 1 - 10 


Mention et al. 




21 


9 - 1050 


~ 3 MeV 


10" 2 - 10" 1 


Bugey 40/15 ratio 


V e 


25 


15 and 40 


3- 8MeV 


> io~ 2 


TOTAL 


V e , Vy, 


48 


10 - 10 a 


3 MeV - 200 GeV 


10~ 4 - 10 a 



TABLE IV: Energies, mixings, and mass splitting sensitivities for each disappearance experiment. 



Unitarity Constraints: The condition for unitarity in the neutrino mass mixing matrix is WU = 1, or 
in component form J2i U a iU^ = 6 a p. In practice, this means that sum of the norms of any single row or 
column in the mixing matrix must equal 1, which bounds the size of any particular element of the matrix. In 
this way, high confidence measurements of mixing angles for v e and with m, 712, and 713 can set bounds 
on the size of the mixings of v & and with 77,4 and 77,5. 

For instance, solar neutrino experiments such as KamLAND [54] measure v e disappearance via the mixing 

sin 2 28 so i = 4|[/ e 2| 2 ^1 — \U e 2\ 2 — \U e sf — \U e 4\ 2 — \U e 5\ 2 ^j. We can extremize this over the mixing \U e 2\ 2 , 
but we find that the limits on |t/ e 4| 2 + |?7 e 5| 2 are not very constraining because the solar mixing angle is 



4 Even though the new reactor flux has been reported to reflect oscillations of a single sterile neutrino, we find that it is not 
consistent with our LSND and MB preferred region for the light sterile neutrino, and we use the reactor data as a constraint. 
On the other hand, it may be fit well with the fifth neutrino of the 3+1+1 scenario, as shown in Sec. IIII Al 
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measured at low confidence and the mixing angle is not maximal. The limits from reactor v e disappearance 
experiments provide stronger limits. 

On the other hand, a similar analysis is effective in constraining 1 1 2 and |f/^5| 2 from the Super- 
Kamiokande data, since the atmospheric mixing angle is measured at high confidence to be maximal: a 
substantial mixing with heavy sterile neutrinos would imply a larger than observed ratio in upward to down- 

2 2 / 2 2 2 \ 

ward going muon neutrino fluxes. This mixing angle is sin 29 atm = 4 |^3| ( 1 — — \U^\ — I \ J , 

and extremizing to find the largest value of |£7^4| 2 + |E/^5| 2 compatible with the measurements gives 

sin 2 20 otm < (l - \U^\ 2 - |C/ M5 | 2 ) 2 ■ (Al) 

Using the global best fit value for atmospheric mixing angle we find the 90% (99%) confidence level con- 
straints: 

lU^l 2 + l^sl 2 < 0.0175 (0.0274). (A2) 

These bounds on 1 1 2 -I- |t^5| 2 are included in the disappearance constraints, and are in practice the 
strongest constraints available. 

Reactor and Gallium Anomalies: The RAA is detailed in (TTJ US] and corresponds to a lower-than- 
expected flux of D e emitted from nuclear reactors. The GA has been reported in and also discussed 
clearly in [T7J and [5SJ: SAGE and Gallex have independently measured a lower-than-expected flux of 
v e from the decay of megacurie sources of 51 Cr and 37 Ar, corresponding to anomalously low rates of the 
reaction v e ± 71 Ga— > 71 Ge + e~ . This is in principle bounded by similar measurements of the rate of 
is e + 12 C —> 12 N g s . +e~ 28; 30 , which are more consistent with expectations [35]. We find that these carbon 
data are compatible with fits to the RAA and GA data. 

We summarize the status of the RAA and GA anomalies as well as the carbon data in Table [V] In Eq. ^ 
we give our fit to the combined data (using correlation information as reported in [IT]). 



Anomaly 


# points 


|C/ e5 | 2 


Xmin 


Gallium 


4 


0.0708 ±0.0317 


1.7 


Carbon 


11 


0.0901 ±0.0874 


8.5 


Reactor 


19 


0.0266 ± 0.0144 


7.2 


TOTAL 


34 


0.0356 ±0.0130 


19.4 



TABLE V: Fits to the reactor, gallium, and carbon anomalies. 



Statistical Methods Employed: To place constraints on neutrino mixing from the various null ap- 
pearance and disappearance experiments we use the raster-scan method described in [56] . For each value 
of Am 2 ^, the \ 2 is minimized with respect to the remaining mixing parameters. For the n remaining 
independent mixing parameters, a A\ 2 test is performed to give an n-dimensional confidence interval at 
each Am| x . In the 3+1 scenario, this corresponds to finding 1-dimensional confidence interval for sin 2 29 at 
each given value of Ato 2 ^. The raster-scan provides a more precise confidence region than a global fit. As 
a result of the sinusoidal dependence of the oscillation probability on Am 2 ^ , the probability distribution for 
the x 2 may deviate from gaussian for large deviations from the true value. This may result both in finding 
an "incorrect" minimum of the \ 2 an d using an incorrect probability distribution function for determining 
the sizes of the confidence intervals. By performing a raster scan, one removes the sinusoidal dependence so 
that the data follow a standard \ 2 distribution. 

Although this a powerful technique for forming exclusions from null experiments, it is less applicable to 
cases in which there is a positive result. This is because the raster-scan does not identify preferred values 
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of the parameter Am^. For this reason we perform a global fit to the LSND and MB data, minimizing the 
X 2 with respect to all parameters. The confidence region is given by a A% 2 for 2 DOF in 3+1 model and 4 
DOF in the 3+1+1 model. 



Appendix B: Oscillation Formalism 



For nonrelativistic neutrinos whose wavepackets travel with same energy E and whose momenta may be 



Taylor expanded as pi — E — m 2 /2E, the probability of oscillation to flavor vn from flavor v a is 

Pu a ^u p = U «iUpiU aj U* 0j exp [i(m? - m))L/2E] 



(Bl) 



This formula is easy to evaluate in the limit of many light mass eigenstates because unitarity simplifies the 
evaluation of the sum. However, the fifth neutrino may either not be accessible or may be produced in a 
reduced phase space. In this case the evaluation of the sum is less straightforward because the assumption 
that production processes for all mass eigenstates are similar may no longer be good. With the definitions 
Uafiij = U^UpiUajUpj and Xij = (to 2 — m 2 )L/4E and taking a phase space suppression factor < a < 1 
on oscillations through 77,5, the oscillation probability is 



= Y,ij U Q/ 3ij exp(2ixij) - 25R (1 - a) V ■ V a p 5l exp(2ixij) + (1 - a)U a/355 



(B2) 



Carrying out some standard simplifications allows us to write Eq. (|B2| as 
P, 



la/3 



[1 - 2(1 - a)|C/ a5 ||C/> 5 |] + (1 - a)\U ab \ 2 \Up b \ 2 - 4£ 



5>i>j 



- 4a Y$=i ^[U a /35j] sin 2 x 5j - 2 J^5>i>j 9pa/3ij] sin 2x tJ - 2a Y$=i $*W a p5j] sin 2x 5j 

(B3) 

where we use the unitarity condition U a iU% i = S a p. In the limit a — > 1 we recover the standard result. 

In all of the experiments of interest we may ignore oscillations due to the mass differences Am 2 2 and 
Am 2 3 , and Eq. (B3) simplifies. The oscillation probability of interest in appearance searches such as LSND 
is found to be 



\U e4 \ 2 \U^\ 2 {o[(l - r) 2 + 4r sin 2 0] + (1 - rf + 4r sin 2 ^! ± /?)} 



(B4) 



where + (— ) is for v (y) oscillations, and the definitions of r and /3 are given in Eq. ^ in the text. For 
disappearance experiments, the relevant formula is 



1-P„ 



sin 2 20 a4 sin 2 x 41 + 2 1 U a5 \ 2 ( 1 - s±l | U a5 \ 2 ) , 



(B5) 



where sin 2 29 a4 = 4|[/ Q 4| 2 (1 — |i7 Q 4| 2 — |/7 a 5| 2 ) and we assume that oscillations through n§ are averaged. 
We emphasize that because the phase space factor a only enters at second order in |C/ a 5| 2 the phase space 
available to 715 has very little impact on the predictions and constraints of disappearance experiments. 

The "zero-distance" effect [33] arises because even in the limits a — > and X41 — > there remains a 
nonzero probability for oscillation. Another consequence of the very heavy state is that for fixed a the 
sum of Eq. ( B3 1 over (3 betrays nonunitarity. This indicates that we have normalized our states incorrectly. 



However, this is cancelled by an inverse change in the production and detection cross-sections in the types 
of experiments considered here |57) . so we may use the formulae as if the probabilities were unitary. 
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